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A-phase origin inn cubic helimagnets (MnSi, FeGe etc) is explained. It is shown that its upper 
bound is a result of the spin-wave instability at H± > ^/8/3A = Ha2 where H± are the magnetic 
field perpendicular to the helix axis k and A the spin-wave gap respectively. The last appears due 
to the spin-wave interaction if one takes into account that the Dzyaloshinskii-Moriya interaction 
acts between different spins. The infra-red divergences (IRD) in the 1/S series for the magnetic 
energy at H± —¥ Ha2 are responsible for the lower A-phase boundary Hai- It is shown that the 
A-phase exists at all T < Tc but if T <C Tc it is very narrow and can not be observed. However in 
the critical region just below Tc its width increasing strongly. Preliminary estimations demonstrate 
semi-quantitative agreement with the existing experimental data. 

The existence of the spin-wave gap is a crucial point of our consideration. In Appendix A we 
' present its derivation in more transparent form than in previous publication. 
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I. INTRODUCTION 



Unusual properties of noncentrosymmetric cubic B20 helimagnets (MnSi, FeGe and related compounds) with 
Dzyaloshinskii-Moriya interaction (DMI) attracted a lot of attention during more than thirty years (see for example^— 
and references therein). Renascence in this field began with a discovery of a quantum phase transition to a disordered 
(partially ordered) state in MnSi at high pressure^—. Then this transition was observed in FeGe also^. 

Apparently complicated behavior of the helix wave- vector k (the helix axis) in magnetic field H is one of the most 
striking phenomenon observed in these compounds. Indeed suppressing the weak cubic anisotropy the field aligns 
the helix axis along itself and gives rise a conical magnetic structure which transforms to the ferromagnetic one at 
critical field He- This simple behavior holds almost in whole region of the (H, T) phase diagram. However just below 
the transition Tq to paramagnetic state the helix vector k suddenly rotates perpendicular to the field. This so-called 
i^i . A-phase discovered by B. Lebech in FeGe^ exists in rather narrow field range Hai < H < Hai < He and above 
Hc2 the helix axis returns to the field direction agai n 1 ! 2 ' 10 : 11 . 

Further small angle neutron scattering experiments revealed unexpected feature. If the neutron beam is along the 
field the small-angle magnetic scattering appears in the A-Phase only and the magnetic Bragg reflections form the 
six-fold structure which dependent weakly on the field orientation relative to the crystal axes. This phenomenon was 
ly-^ observed in MnS&, FeGeH and several B20 compounds 1 . It could be described as superposition of three helices with 
£f~) \ wave vectors ki + k2 + k3 = where |k,| = k which are perpendicular to the field H. 

Meanwhile it was claimed that this structure is a new skyrmion lattice state 1 ' 2 . However the higher-order reflections 
which has to be in this case were not observed. Hence a nature of this six-fold structure inside the-A has not been 
understood yet. 

Moreover the very existence of the A-phase with perpendicular k orientation has to be explained. Indeed in zero 
field the multi-domain state is realized with k-axes are along (111) or (001) depending on a sign of very weak cubic 
anisotropy 13 . Then with the field increasing the vectors k rotate to the field and conical helix structure develops 
with the cone angle determined by sin a = —H/H o 14 ' * 15 - This behavior is the same as in antiferromagnets above 
the spin-flop transition where we have not any hint to the A-phase state. It should be noted also that the classical 
magnetic energy depends on H|| the field component along the helix vector k only and Hj_ dependence appears as a 
quantum phenomenon onlyis. 

Behavior of the cubic helimagnets in the field is related with more important multiferroic problem. Indeed in 
RMnOs materials multiferroic properties are connected directly to spin helices mediated by DMIi£~— and a lot of 
transitions in magnetic field were observed (seei^ and references therein). However we have not now any satisfied 
explanation of these transitions. Hence understanding of the helices behavior in magnetic field may be considered as 
an urgent problem. 

In this paper we explain the A-phase origin developing theory of the helix behavior in the field which takes into 
account the interaction between spin- waves. This interaction gives rise the infra-red divergences (IRD) in perturbation 
1/5 series for the field-depending part of the magnetic energy at H\ — iJ| 2 — »• 0, where is the field component 
perpendicular to the helix vector k and Ha2 is the upper A-phase boundary respectively. We demonstrate also that 




FIG. 1: Two parts of the magnetic energy En and E± with the the helix vector k along and perpendicular to the field 
respectively. Between Hai and Hai in the A-phase region present approach is applicable in the multi-domain case only. 



Ha2 — where A is the spin- wave gap. In the Hartree-Fock approximation it is given by^s 

where S is the unit-cell spin. 

We demonstrate that the magnetic energy is a sum of two parts Em = Eu + Ej_ differently depending on Hn and 
H± as shown in Fig.l. This anisotropy is a consequence of an unusual form of the spin- wave energy at Hn <C He 
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^{AkY{q\ + 3g 4 /8fc 2 ) + A 2 - 3HI/8, q < k; 



MVl 2 + k 2 , q > k, 



(2) 



where A is the spin-wave stiffness at q 3> k. This q anisotropy is a result of the DMI which mixes spin-waves with q 
and q ± k~i2i. These umklapps produce satellite spin- wave structure with q ± nk where n = 1, 2, ... observed in 2 ^. 

It should be noted that in 21 - was claimed that in B20 helimagnets the spin-waves are the gapless Goldstones due to 
translation invariance along the helix vector k. Really this statement is correct in the linear theory only and is broken 
in higher approximations if one takes into account that the DMI acts between different spin o 14 ' 23 " — . The essence of 
the problem is following. Conventional macroscopic expression for the DMI energy is given by Z?(M(r) • [V x M(r)]) 
is ill-defined as both M operators act in the same r point. It is unimportant in the cases of the classical groundOstatc 
energy and the linear spin-wave theory but becomes crucial if one is interested in the spin-wave interaction. Indeed 
if both M operators act in the same point they non-commute and we have not a gap. Otherwise they commute 
and the gap appears in the Hartree-Fock approximatio n i 23 (see also line below Eq.(15) and Appendix A). In this 
case the DMI feels the lattice structure the above mentioned translation invariance is broken. The sum in Eq.(l) 
for A 2 is saturated at q *~ 1/a where the macroscopic approach is not applicable. Moreover in 2 ^ was shown that 
the magneto-elastic interaction mixes k = magnons with the 2k phonons. As a result in the second order on this 
interaction a negative contribution to A 2 appears and we have 

A 2 = A 2 , F + A 2 f£ , (3) 

It was suggeste d 24 ! 25 that the quantum phase transition at pressure observed in 5 .^ is a result of a competition 
between these two terms and holds when A 2 = 0. 

From Eq.(2) follows also that uniform spin- waves are unstable if 



H ± > H A2 = AV8/3 



(4) 
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and the perpendicular state is impossible. Hence Ha2 is the upper A-phase bound. Below in agreement with 
experiment we assume that H A2 < H c and demonstrate that the lower bound Hai is a result of the infra-red 
divergences (IRD) which appear in the 1/5 perturbation expansion for the magnetic energy if H± —> Ha2- At low 
T this boundary is very close to Ha2 and the A-phase can not be seen. However just below the transition Hai 
decreases strongly due to critical slowing down (decreasing of the spin- wave stiffness A as T — > Tc) and the A-phase 
becomes visible. Qualitative agreement with experimental data of Ref.— was demonstrated. For detailed comparison 
with experiment one has to have more precise experimental data for the T dependence of Hai,2 in the critical region. 

II. SPIN- WAVE INTERACTION 

At the beginning we have to summarize briefly principal theoretical results which will be explored below. We 
use the Bak- Jensen models adding the Zeeman energy and omitting weak cubic anisotropy which is unimportant. 
Corresponding Hamiltonian is given by 

H = J q( s q ' s -q)/2 + *A,(q • [S q x S_ q ]) + Vn (H ■ S„)}, (5) 

where the first second and third terms are the ferromagnetic exchange interaction, DMI and the Zeeman energy 
respectively. In the case of the helical structure for S q we havei 4 - 

S q = SZc + SlA + Sf A* 

. (6) 
Sq = sin a + 5 q cos a, 5 q = S^_ k cos a — 5* q _ k sin a + iS2^_ k , 5 q = S^ +k cos a — Sq +k sin a — iS^ +k) 

where A = (a — ib)/2, unit vectors d, b, c form the right-handed orthogonal frame and the spin operators in the 
Dyson-Maleyev representation are given by 

5 q = 7V 1/2 S<5 q , - (a+a) q ; S* = -iy/S/2[a* - at q - (a+a 2 ) q /2S] ; S* = y/S/2[a^ + a± q - (a + a 2 ) q /2S], (7) 

where a q and a q are conventional Bose operators. 

In the classical approximation we have k = SD^c/A where A = S(Jq — J k )/fc 2 is the spin- wave stiffness at q 3> k, 
sin a = —Hu/Hc where H\\ is the field component along k and He — Afc 213 ' 14 . Corresponding part of the magnetic 
energy have the form E{H\\) = —H 2 /2Hc- 

We are interested below low-field region H < He and for simplicity put a = 0. Conventional spin- wave Hamiltonian 
is given byi^ 

H 2 = ^[S q a+a q + S q (a q a_ q + at q a+)/2], (8) 

where at Hu <C He we have 1 ^ 

= S(M 0;k -M q , u ) + J B q ~ A{q 2 + k 2 /2),B cl = S(M q , k - J q )/2 ~ Ak 2 /2, Af q , k = (J q+k + J q _ k )/2 + 2 J D q (k-c), (9) 

where approximate equalities hold at q <C 1/a and a is the lattice spacing. The equilibrium condition is given by 
£>o(k • c) = Ak 2 /S. 

The second term in Eq.(5) mixes excitations with q and q ± k and gives rise the q anisotropy in Eq.(2) at q < fc 14 ' 21 . 
For following we has to consider the spin- wave interaction in more details in comparison witbM. From Eqs. (5-7) at 
small a for the interaction energy we have Vi = V4, + Ve and 

Vi = -(l/2)^(A/i_3,k-Afi, k )a+at 1+2+3 a 2 a3 + (l/4)^(Ji-Afi. k )(ai+a^ 1 )(a + a 2 )-i, (10) 
V 6 = -(1/165) ^(Ji - Af 1 , k )(a+a 2 ) 1 ( a +a 2 )_ 1 ~ (Ak 2 /16S 2 ) ^(a+a 2 ) 1 (a+a 2 )_ 1) (11) 

where 1, 2, 3 label corresponding momenta. The first term in Eq.(10) is a generalization of the Dyson interaction for 
ferromagnets 2 ^. We do not consider here the Ve interaction as it gives further terms in 1/S expansion. 
In the Hartree-Fock (HF) approximation Eq.(10) gives rise non-Hcrmitian spin- wave Hamiltonian 

Hsw = ^[£ q a+a q + (B+a q a_ q + S q al q a+)/2], (12) 
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where Eq = Eq + E, -B^" 1 = B q + n^ + ) and corresponding diagram is shown in Fig. 2a. Corresponding expressions at 
q = are given by 

E = (1/2) Y^[{Ji ~ Af llk )(m + /i) + (Jo - M ok )(m + /i/2)]; (13) 
n = ]T[(M , k - Mi, k ) + i(J - Mo, k )]/i; 
n+ = (1/2) ^[(Ji - Afi lk )(m + /i) + 2(J - Af 0jfc )ni]. 

where 

»i =< a+ai >= [(2iVi + l)^i - eJ/26!, /i = -^(2^ + l)/2 £l (14) 

where iVi is the Plank function. From these expressions using definitions (9) for the spin-wave gap we obtain (for 
discussion see Appendix A) 

A% F = 2£ E - B Q (U + H+) = J2 (15) 

where Eq = Bq = Ak 2 /2. We have Aj^ F ^ due to —1/2 term in expression for n%- Forgetting that the DMI acts 
always between different spins we must replace m — > n% + 1/2 and A 2 HF = 02£. It should be noted also that the sum 
in Eq.(13) is saturated at q ~ 1/a and may be considered as T-independent at T < To- 
One has to note also that Fig. 2a diagram is of order of S° and A 2 ~ S whereas the main part of e 2 ~ S 2 . 

III. PERPENDICULAR FIELD (H ± k) 

From Eqs.(5-7) for interaction with the perpendicular field at a -C 1 we have^ 

V x = (H • A)[y/S/2(*-* - a+) - <_k«q] + h - c - = V l + Vum, (16) 

where the terms linear in the operators lead to the spin-wave Bose condensation at q = ±k. Other terms mix 
excitations with q and q ± k 

in Appendix B we ob 
their solution we have 2 ^ 

a k = VW2(H • A)/Ak 2 ; a_ k = -^/S/2(H ■ A*)/Ak 2 . (17) 

and the magnetic energy is given by»22. 

SH 2 SHl 

Em °--2Hc-~4Hc-- (18) 

where the factor 1/4 is very transparent: 1/4 = 1/2 < cos 2 ip > where ip is the angle between the helical spin and 
H^. According to this equation the helix axis k has to be along the field as was observed in all H,T region of the 
phase diagram except small A-phase pocket just below Tg ^ 10 ' 11 . 

Taking into account the spin- wave BC we get from V4 additional contribution to the spin- wave Hamiltonian (12) 

H BC = -(Hl/lQAk 2 ) J2(a+aq + a+at q + 2a q a_ q ). (19) 

As a result we obtained —3Hj_/8 term in Eq.(l) instead of —H\/2 inM This part of the spin- wave Hamiltonian has 
been taken into account in Appendix B for the Green function evaluating and produce non-symmetric H\ terms in 
Fq and F+ Green functions (see below). 

These umklapps gives rise to two infinite systems of conjugated linear equations for the Green functions. At 
H± <C He they can be truncated and we have two systems of the linear equations considered in Appendix. Their 
solution may be divided on two parts: direct and umklapp Green functions. For the first one we have 

= iu + Eq + V-7H 2 JlQAk 2 _ Bg + U Bq + B+ - Hj /8Ak 2 

Gq{u) - Gq{-U) -j— 2 , Fq - -j-j-^ , F q - . (20) 



in Appendix B we obtain a system of equations for c&q and at^ considering them as classical variables. From 
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FIG. 2: Diagrammatic expansion for the magnetic energy. Lines are Green functions G =<—, F and F + =< >. At top 

is the sum of the HF loop and the BC part Eq.(19) used in the diagrams. The IRD first time appears in Fig.2d diagram. 

The umklapp functions are given by 

G+ = a - = ' Ft = ' F - = '^rkr G -^ a+ = ' F - = ' F+ = '^T^r (21) 

where e q is given by Eq.(2)and in numerators we neglected some terms proportional to q\ and q\ which are unim- 
portant to us (see 14 ). Expressions (14) for n q and / q are direct consequences of these equations. 

IV. MAGNETIC ENERGY 

We demonstrate below that the A-phase lower bound is a result of the infrared divergence (IRD) at Ha2 — H±_ — > 
which appear in the 1/S perturbation series for the interaction energy. Using Eqs.(12, 18,19) for the total magnetic 
energy we have 

E = E M o+<Hsw + H bc > + <V-Vhf>, (22) 
where the first term us given by Eq.(18). For the second term we have very transparent expression 

Esw = $> q iV q + (e q -.E q )/2], (23) 

where off-diagonal terms in Eq.(16) are responsible for correct —3Hj_/8 field dependence of e q . This expression 
includes the energy of thermally excited magnons and zero-point motion with energy given by Eq.(2). The last term 
of Eq.(22) does not contain the HF part Fig. 2b as it was used in Eq.(12) ant the 1/5* expansion begins with Fig. 2b 
diagram. 

In the principal 1/S order left and right vertexes in Fig. 2 IRD diagrams may be expressed using Eqs.(13) 

r = [e - (n + n+)/2] a + flq + (i/2)(n<z+ + n+ a _ q )(a q + oi q ), (24) 

where according to Eqs.(19,B5)we have to add to E,I[( + ) the BC parts (fir st line in Fig. 2). . As a result we get 

T = (A 2 /Afc 2 )[i? ia +a q - (R 2 a+ + i? 3 a- 9 )(a q + <)/2)] (25) 

where Ri = 1 + H\jm\ 2 ~ 4/3, R 2 = 1 + H\/QH 2 A2 ~ 7/6, and R 3 = H\jZH\ 2 ~ 1/3. Approximate equalities 
hold et H± — > Ha2 and will be used below. 

Using Eq.(25) and taking into account direct (20) and umklapp (21) Green functions for the first IRD diagram 
Fig. 2d we obtain 

_ H\ 2 T v l + gj/32Afc 2 
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where we neglected in the numerator terms bilinear in (Hj_, A 2 ) and in following will omit the Hj_ term. 
We are interested in T Ak 2 region and the IRD displaces the oj n = term only and we obtain 

T(ak) 3 ^/3/8H A2 ( Hi \ _ H\ . . 

E d = 75T7TT2\2 775 EFT + 1 - ^772-' ( 27 ) 



where a is the lattice constant, and instead of Eq.(18) we have 

SHI SH\ ( AH c E d 



e --u&-3&\?sb£) ™ 

Neglecting higher order IRD terms (see below) we obtain that the A-phase is energetically profitable if the expression 
in the brackets becomes larger than two. Hence at the S-phase lower boundary we have 

E d = -4SH 2 A2 /H C (29) 

and using Ea.(27)we get 

Hax = H A2 [1 - 2^/3/8(ka) 3 TH c /9TTS(Ak 2 ) 2 } 1 ^ (30) 

Using well known low-T parameters for MnSi: T c ~ 29.5K, H c ^ Ak 2 ~ 0.6T ~ 0.8K, a ~ 0.46nm, k ~ 0.3871m" 1 
and S = 1.6(see for example 1 ^ and references therein) we obtain H A i/H A z = (1 — O.OOSST/Tc) 1 / 2 . Hence at 
T <C Tc the A-phase exists but it is so narrow that can be hardly observable. Moreover at T = Tc the ratio 
r = H A1 /H A2 = 0.9975. 

We demonstrate now that there is strong decreasing of the ratio r = Hai/Ha2 just below Tc- It is connected 
mainly with the spin-waves critical slowing down. Unfortunately there are not any theoretical predictions for it and 
our discussion is restricted by analysis of Eq.(30) where all parameters may be measured independently. 

First of all the low-T condition He — Ak 2 is violated. Indeed Hc(Tc) ~ i?c(0)/3 2 -, k is almost T independent^ 
whereas the spin-wave stiffness A is renormalized strongly. Indeed it is equal to 0.52meUnm 2 and 0.24meVnm 2 at 
T = 5K and 26K respectively^^ while H c remains almost unchanged 2 . One may await that near Tc A = Aqt z 
where r = (T c -T)/T c . 

The unit cell spin S — 5qt - 2 228 ' 29 . According to^S this scaling behavior begins at T w 25K (r 0.14) and 
So ~ 1.2. The upper A-phase bound is almost r independent 2 -^ and from Eqs.(l) and (4) we have A ~ S 1 / 2 and 
A = Aot°- u . Putting A a = 0.2AmeVnm 2 we obtain 

H A1 = H A2 (l - 0.033Vt 44 ) 1/2 , (31) 

where h = H c (T)/H c {0). 

The most detailed available data for the A-phase in MnSi are shown in Fig.l of Refj£ and we compare our results 
with them. i. T = 25K, t = 0.14, h ~ 1 r = H A1 /H A2 = 0.96 and A-phase was not seen. ii. T ~ 28.5K, r = 
0.03 h ~ 0.8; r — 0.94 whereas the observed ratio is 0.6 — 0.7. Hence the theory explains the A-phase phenomenon at 
least qualitatively. It has to mention also that all parameters including values of r used above are known with very 
low accuracy. For example at r = 0.01 we have r = 0.89 and replacing in Eq.(31) 0.033 — > 0.1 we get r = 0.63 We 
mention also that according to 2 -^ 

H A2 ~ 0.22T; and A ~ 0.13T = 15fieV. (32) 



V. 1/5 CORRECTIONS 

The next IRD term is represented by Fig.2e diagram and in general form is given by 

E e =< ru 4 r > . (33) 

Interaction V4 given by Eq.(10) consists of two parts. The first one gives the central vertex in Fig.2e diagram which 
disappears at zero momenta and leads to rather weak IRD singularity which may be neglected. The second remains 
constant and gives main correction to the above results. 

Using Eqs.(2),(20) and (A3) for the most singular e2 correction we obtain 

3Ak 2 fA 2 R 1 \ 2 x ^B 1 (F+ -F ± )B 2 15E d Ak 2 H 2 ±x ^T 
Ee --—{-Ak^J ^ Z^l > 32^ ( } 
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where Zi = cl>, 2 + e 2 and I = 1,2. This expression is proportional to {H\ 2 — Hj_) 2 Using Eq.(26) we obtain 
E(r/e 4 ) = SE d /H\ 2 and 

_ 15£ 2 Afc 2 g 2 

^ e - 45i?l 2 (35) 

Considering this expression as a small correction to Ed in Eq.(29) we must replace 

£ d ^£ rf (l + 15^fc 2 /16ff c ). (36) 

At low T this correction of order of unity and one has to examine all 2/5 series. However near Tc where A — > Aot 
it may be neglected. In any case deep into the A-phase and near its upper boundary the full series examining is 
unavoidable if one assume conventional helical structure (c£). 



VI. CONCLUSIONS 



A-phase origin in cubic helimagnets (MnSi, FeGe etc) is examined. It is shown that its upper bound is a result of 
the spin- wave instability at H±_ > ^/8/3A = Ha2- 

The infra-red divergences (IRD) in the 1/5 series for the magnetic energy at H± — > Hai are responsible for the 
lower A-phase boundary Hai- It is shown that these IRD are anomalously strong due to softens of the spin- wave 
spectrum at q < k. As a result the low- momenta fluctuations behave as in 2D systems^. 

We demonstrate that the A-phase exists at all T < Tc but if T <C Tc it is very narrow and can not be observed. 
However in the critical region just below Tc its width increasing strongly. Preliminary estimations show semi- 
quantitative agreement with the existing experimental data. However the problem demands further theoretical and 
experimental studies. We can formulate following unresolved problems, i. Examination of the full 1/5 expansion at 
low T as the second correction to the Hai boundary is not small. In this respect we wish to point out to the hint to 
the A-phase observed in MnSi at T = 10K (Fig. 3b in Ref. 15 ). ii. Better understanding the temperature dependence 
of the spin-wave energy at q < k. iii. More precise measurements of the A-phase boundaries then it was done in 
previous studies. In this respect we wish to note that there are the first order transitions at both a-phase boundaries 
which has to be accompanied by the specific heat jumps. The first time the specific heat anomalies in this region were 
observed in^si but more precise measurements would be important. 

In this discussion we avoided the nature of the A-phase itself. There are two possibilities: exotic "skyrmion" state 
proposed in 2 - or the tree domain structure. However at present we have not any real theory of the skyrmion state as 
well as an explanation of the three domain state. 
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Appendix A 

As very existence of the gap is crucial for us we repeat here some details of corresponding calculations which were 
presented inM in a rather cumbersome form. In Eq.(15) terms proportional to Jo — M q .k cancel. For remained parts 
of S, n(+) using Eqs.(9,13,14) we obtain 

E = n+; n = E - ^(M llk - JO/4 (Al) 

£(M llk - JO/4 = (1/2) ££> q (k • c) (A2) 

as ^2 Ji ; k = E J q = due to condition that the exchange interaction acts between different spins. Using equality 
L>o(k • cj = Ak 2 /S we obtain Eq.(l). 

It may be shown that expression (A2) is much larger than other terms in Eqs.(13) as they have additional small 
factors which vanish if k — > 0. Hence we may put 



n = -2A 2 /Ak 2 , E = n+ = 0. 



(A3) 
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From Eq.(A2) follows that A 2 7^ if Dq 7^ 0. We demonstrate now that this condition is fulfilled. for the 
DMI was suggested expression 

r,R'(V - V')[Sr x Sr/] = i^D q (q • [S q x S_ q ]), (A4) 

where Drr' = Dr».r, the condition R^R' holds and ^2 q-D q =0. 
We present now simple example where the ^2 Dq 7^ was evaluated. 

In2S a simple model for the DMI in cubic lattice was proposed where the DMI acts between neighboring spins and 
the DMI vector dj t q — —iK svaqj, j — x,y,z.. For the form-factors we have (si&qj)/qj and sin qj / qj = ir/2. 

It should be pointed again the gap appears if the DMI acts between different spin and the translation invariance 
along k is broken. Otherwise this invariance is restored and the gap disappears due to non-commutativity of two spin 
operators in the single lattice point22. 



Appendix B 

At H± <C He truncated equations for the Bose condensed magnons are given by (cfj^) 

E a+ + B+a -ha+- fa\ = 0, (Bl) 



Ba£ + E a 


- /ia_k - /a k = 


0, 


-f4 + 


Eia£ + -Ba_k = 




-fao + 


Ba^_ + i?ia_ k = 




—hciQ + 


Eiat k + Bak = 


-y/S/2h, 


—ha + 


£?a+ k + Si ak = 





and we obtain Eqs.(17). 

We have two sets of the Green functions determined as follow 

Gq = - <Taq,a+ >;F+ = - <Tat ci ,a+ >;G ± = - <Taq ±k ,a+ >;F+ = - <Tat clTk ,a+ >, (B2) 
G q = - < Ta q ,a q >;F q = - < Ta- q ,a q >;G± = - < ra q±k ,a q >;F± = - < Tat,^, a q >, (B3) 

They are solutions of two conjugated systems of linear equations. W e consider here the second set only. Corresponding 
truncated system is given by 

(iw + E)G + B+F - hG+ - fG- = -1, (B4) 

—BG + {iu: - E)F + hF+ + fF- = 0, 

—fG + (iuj + £i)G+ + BF+ = 0, 

fF - BG+ + (ico - E 1 )F+ = 0, 
-hG + (ico + E 1 )G- + BF- = 0, 

hF- - BG- + {iu: - £1) = 0, 

where h = f* = (A • H), E± ~ S q ±k if 5 fc and according to Eqs.(12,13) and the BC contribution (19) we have 

E = Eq + S - Hl/16Ak 2 ; B = B + U- Hl/8Ak 2 ; B+ = B + U+ - H\/AAk 2 , (B5) 

Solution of these equations as well as the conjugated one are given by Eqs.(20) and (21). 
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